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The Hubbard model on the Kagome lattice is investigated in a metallic phase at half-filling. 
By introducing anisotropic electron hopping on the lattice, we control geometrical frustration and 
clarify how the lattice geometry affects physical properties. By means of the fluctuation exchange 
(FLEX) approximation, we calculate the spin and charge susceptibilities, the one-particle spectral 
function, the quasi-particle renormalization factor, and the Fermi velocity, ft is found that geomet- 
rical frustration of the Kagome lattice suppresses the instability to various ordered states through 
the strong reduction of the wavevector dependence of susceptibilities, thereby stabilizing the for- 
mation of quasi-particles due to the almost isotropic spin fluctuations in the Brillouin zone. These 
characteristic properties are discussed in connection with the effects of geometrical frustration in 
the strong coupling regime. 

PACS numbers: 71.10.-w, 71.27. +a 



I. INTRODUCTION 

Geometrically frustrated metallic systems have at- 
tracted much interest since the discovery of the 
heavy fermion behavior in the transition metal oxide, 
LiV^O^^. This compound has the spinel structure and 
the band calculations show that the conduction bands 
are essentially composed of vanadium 3d orbitals, which 
are well separate from the oxygen 2p bands. There- 
fore when low-energy properties are concerned, it is suf- 
ficient to consider only the vanadium 3d orbitals 3-4-5-6 . 
The vanadium sublattice in the spinel structure consti- 
tutes a network of corner-sharing tetrahedra, i.e., the 
pyrochlore lattice, and this is a typical example of ge- 
ometrically frustrated lattices in three dimensions (3D). 
Various works have revealed unusual low-energy proper- 
ties of spin systems on the pyrochlore lattice, including 
the presence of thermodynamic degeneracy of the ground 
states£*&2ii2iiiii2ii£. Another metallic pyrochlore com- 
pound, Y(Sc)Mn2, also exhibits the heavy fermion be- 
havior where unusual feature has been observed in the 
dynamical susceptibility 14-15 . In both compounds, no 
long-range order has been observed and the specific heat 
coefficient is strongly enhanced at low temperatures, sim- 
ilar to the lanthanide or actinide heavy- fermion systems. 

The heavy fermion behavior in the lanthanide or ac- 
tinide systems is, as well known, attributed to the Kondo 
effect and the presence of localized /-orbital is essential 
to this mechanism 1 -. However, in these transition-metal 
heavy fermion systems, 3d electrons are much more mo- 
bile than /-electrons and the presence of local magnetic 
moment has not been detected so far. It is highly non- 
trivial whether the enormous mass enhancement is also 
attributed to the Kondo effect, and there may exist an- 
other mechanism leading to the heavy fermion-like be- 
havior. This indeed gives rise to a number of theoretical 
proposals on the mechanism of the formation of heavy 



quasi-particles. Among them, it has been claimed that 
geometrical frustration plays an important role even in 
the metallic phase when the electron interactions are suf- 
ficiently largo 17 i 18 i 19 i 20 . 21 i 22 i 23 i 24 i 25 i 26 i 27 i 28 . 

The two-dimensional (2D) Kagome lattice is another 
geometrically frustrated system, which may be regarded 
as a 2D analog of the pyrochlore lattice. Antiferromag- 
netic spin systems on this lattice have been intensively 
studied and many unusual properties have been found. 
For example, in the 5=1/2 Heisenberg antiferromag- 
net, there exists a finite energy gap between the sin- 
glet ground state and triplet excitations, and a thermo- 
dynamic number of singlet excitations exist within the 
singlet-triplet gap due to strong frustration22i22*2ii22i2£. 
This level scheme of excitations is in common with that 
for the spin systems on the 3D pyrochlore lattice, and it 
may be quite generic in geometrically frustrated systems. 
If electrons become itinerant, electron motion will be 
coupled to both spin-triplet and spin-singlet excitations. 
Therefore, quasi-particles in this case may be renormal- 
ized in a different manner than other more conventional 
cases without frustration. To study this issue, we will 
investigate the correlation effects in the 2D Kagome lat- 
tice. The Kagome lattice is simpler than the pyrochlore 
lattice due to the low dimensionality, but otherwise its 
magnetic properties in the insulating phase are in com- 
mon with those in the pyrochlore lattice. Therefore, the 
study on a metallic Kagome system is a good starting 
point to investigate how geometrical frustration affects 
physical properties in a metallic phase, and we expect 
that many of the results will hold in the other frustrated 
systems including the pyrochlore system. 

In this paper, we will investigate the effects of geomet- 
rical frustration of the Kagome lattice on physical proper- 
ties. In particular, we focus on the possibility of magnetic 
instability and its relationship to the nature of quasi- 
particles in a metallic phase. We employ the fluctuation 
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exchange (FLEX) approximation 36 i 37 i 38 i 39 i 40 i 41 for elec- 
tron correlations and calculate the spin and charge sus- 
ceptibilities, and the one-particle spectral function. We 
will show that geometrical frustration of the Kagome lat- 
tice indeed suppresses the instability to various ordered 
states, and stabilizes the formation of quasi-particlcs. 

This paper is organized as follows. In the next section, 
the model and the method are described briefly, and we 
show the obtained results in Sec. III. Brief summary is 
given in Sec. IV. 



II. MODEL AND METHOD 

The original Kagome lattice is schematically shown in 
Fig- d ( a )j which is given by a corner-sharing 2D net- 
work of triangles. This is topologically equivalent to the 
decorated square lattice with specific diagonal bonds, as 
shown in Fig. d (b) , where the unit cell contains three 
sites. In the following discussions we will deal with the 
latter lattice, which makes the analysis simpler, since its 
Brillouin zone (BZ) becomes from hexagon (Fig. Q](c)) 
to square, —it /a < k x ,k y < n/a (Fig. [T] (d)), where 
a is the lattice constant. In order to control geometri- 
cal frustration on the Kagome lattice, we here introduce 
anisotropic hopping, t and t', between nearest neighbor 
sites as shown in Fig. d (b). When t'/t = 1.0, the sys- 
tem is equivalent to the original Kagome lattice. The 
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FIG. 1: (a) Original Kagome lattice, and (b) decorated square 
lattice which is topologically equivalent to (a). Lattice points 
labeled by m — 1, 2, and 3 represent relative positions in each 
unit cell, t and t' are hopping integrals, (c) and (d) BZ of the 
original Kagome and decorated square lattices, respectively. 
A-D represent the corresponding wavevectors in both lattices. 
Hexagon shown by dash-dotted line in (c) is the 1st BZ of the 
Kagome lattice, which is equivalent to the diamond shown by 
dotted line. Square shown by dotted line in (d) represents the 
1st BZ corresponding to the diamond in (c). 



advantage of introducing this anisotropy is that it prop- 
erly interpolates the non- frustrated lattice (small t'/t) to 
the fully-frustrated lattice (t'/t = 1) in the strong cou- 
pling limit U/t — > co at half filling, where the system is 
reduced to the antiferromagnetic Heisenberg model. We 
expect that the effects due to the frustrated lattice ge- 
ometry should show up even at weak or intermediate val- 
ues of U/t. Hereafter we call the decorated square lattice 
with t'/t — 1.0 (t'/t =^ 1.0) isotropic Kagome (anisotropic 
Kagome) lattice and the lattice constant a is taken to be 
unity, a = 1. 

Then we consider the Hubbard model on the 
isotropic/anisotropic Kagome lattice of Fig. [2(b). The 
Hamiltonian is given by 



H 



£ k C kmcr C km'cr A* 2_v C kmcr C kira<x 
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+ ^ £ C imT Cim T C iml C 



iml j 
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where Ck mff ( c k mcr ) represents the annihilation (creation) 
operator of an electron at wave vector k with site index 
to and spin a. Note that the lattice points are labeled by 
the position of each unit cell % together with the relative 
position m(= 1, 2, 3) in the cell, there exist three bands. 
Here n is the chemical potential and U is the on-site 
Coulomb repulsion. The kinetic energy e^ m is given in 
the matrix form. 
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In the following discussions, one hopping integral is 
fixed and taken as units of energy, t = 1. In the present 
study, we focus on the half-filling case, in which geomet- 
rical frustration in the strong coupling regime is most 
prominent and discuss how the system changes its charac- 
teristic properties from the weak coupling to the strong- 
coupling regime. 

In order to study electron correlations in the model, 
we employ the FLEX approximation, which is a self- 
consistent perturbation method with respect to the 
Coulomb interaction U. The FLEX approximation is a 
conserving approximation based on the idea of Baym and 
Kadanoff 34 ^, and has been successfully used to describe 
electron correlations in the high-T c cuprates and other 
correlated electron systems^^LSS^iiMi. In our three- 
band system, it is convenient to represent the Green func- 
tion G and the effective interaction V in the 3x3 matrix 
form corresponding to 3 sites in the unit cell. As far 
as there is no spin order, the Green function and other 
quantities are diagonal in spin space, and we drop the 
spin index. The self energy within the FLEX approxi- 
mation is written as 



S(fc) = ^£*K(-g)G(fc + g ), 



(2) 
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where T is the temperature, k = (k, iu) n ) and q = (q, iv{), 
with ui n = (2n - 1)ttT and i/j = 2ZttT, and N' = N/3 is 
the number of total unit cells, with N being the number 
of total sites. The effective interaction is given by 



V(q) = Ui + U 2 x{q) + \u 2 x(q)[(l-Ux{q)) 
+ lu 2 x(q)[(i + Ux(q)) 



x(q) 



N 



T J2G(k + q)G(k), 



-I 
(4) 
(5) 



where / is the unit matrix. The Dyson equation for the 
rcnormalized Green function reads 

G(fc)- 1 - g{k)~ x - t(k), (6) 

where g(k) is the bare Green function defined as 



(iu> n + fi)I - e k 



(7) 



By numerically iterating the procedures of Eqs. I©-®, 
we obtain the renormalized Green function. 

In order to investigate the spin and charge response of 
the system, we introduce the generalized susceptibility 
X'(q) of which element is defined by 

Xmim2,m3m4 (q) 

dre WiT <p mim2 (q,r)pt l3 „ l4 (q,0)), (8) 

'o 

where /w(q,r) = [l/VW) £\ e iq - Ri cf m (T)c im / (r) is 
the generalized polarization at imaginary time r. 
represents the position of the unit cell i. It should be 
noted that generalized susceptibility x' is m the 9x9 
matrix form and is diagonal for the spin sector. Within 
the FLEX approximation, the generalized susceptibility 
is obtained as 



X 



mi?7i2 ,77137714 



(q) 



t [k + q)G m3mi (k) 



-ik-(r mi -r m3 ) i(k+q)-(r„ 



(9) 



where r m (m = 1-3) is a relative lattice position in each 
unit cell (Fig. H^b)). The spin and charge susceptibil- 
ities, x^ an d X j are the response of spin-triplet and 
spin-singlet polarizations, respectively, and their matrix 
element is given by 



y (s,c) (a) 
X'(q)(lTUx'(q) 



mim,2 ,77137714 



(10) 



where the =p sign should read — for the spin and + for 
the charge susceptibilities, respectively. The element of 
U is 



mim2 ,7713777-4 



US rni ^ m2 <5r7l3 ,7714 ^771 1 ,7713 * 



(11) 



III. RESULTS 

We numerically iterate the above-mentioned proce- 
dure of Eqs. until the calculated G(k)'s converge 
within desired accuracy. The summations are efficiently 
carried out by using the fast Fourier transform (FFT) 
with N' = 64 2 points in the k-summation and 2048 Mat- 
subara frequencies in the w ra -summation . 



A. Free Electron Properties 

We start with the non-interacting case, U/t = 0. By 
using an orthogonal transformation, the kinetic term of 
the Hamiltonian is diagonalized at each k. The result for 
the isotropic case t'/t = 1.0 is given by 



1±A/1 



21 



tto (T)-(?)-m 

(12) 



where + and — sign correspond to the lowest band, -Eki, 
and the middle band, , respectively. While the largest 
eigenvalue forms a flat band over the whole BZ, the lower 
two dispersive bands are symmetric with respect to ui = 
— 1.0, and touch each other at k = ±(27r/3, — 27r/3) with 
linear dispersions. Let us define the density of states 
(DOS) with including the chemical potential, as D(ui) = 
(1/N) J2k q l/( tJ + A* — Eka) and they are shown in Fig. 
121(a) for several t' . For the isotropic case t'/t = 1.0, a 
(5-function peak appears at u> ~ 1.5 {u> + fi = 2) due to 
the flat band E\&. In the anisotropic case, the dispersion 
of the three bands is modified, and the main change is 
that the highest band is now dispersive, leading to the 
broadening of the 5-function part of the DOS. However, 
the qualitative change in the DOS around the Fermi level, 
u> = Q, is quite small. Shown in Fig. |3 (b) is the Fermi 
surface for the corresponding cases. When t'/t = 1.0, 
the shape of the Fermi surface is most isotropic. As t' 
decreases, the Fermi surface is elongated along the k y = 
k x direction, and becomes more anisotropic. 
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FIG. 2: (a) DOS and (b) Fermi surfaces of the non-interacting 
system (U/t = 0) at half filling for various hopping amplitude; 
t'/t = 1.0 (solid line), 0.8 (dashed line) and 0.6 (dotted line). 




FIG. 3: (a) Largest eigenvalues of the static susceptibility 
X'(q,0) at half-filling for U/t = 0. (b) Contour plot of (a). 
The circles indicate the location of the peak values. The tem- 
perature is T/t — 0.1. 



Once the free electron Hamiltonian is diagonalized, it 
is easy to calculate the susceptibility using Eq. JUJ). Since 
the unit cell contains three atoms, there are 9 indepen- 
dent modes of polarization concerning the site degrees 
of freedom in the unit cell, and each of them has both 
spin-singlet and spin-triplet branches. Correspondingly, 
the susceptibility is a 9 x 9 matrix at each q, and its 
largest eigenvalue is the dominant response. As we are 
considering the noninteracting case (U/t = 0) at the 
moment, the response of spin-singlet polarization is al- 
ways degenerate with that of spin-triplet polarization. 
Of course, they start to differ, upon switching on the 
electron correlations, which will be discussed in a later 
part. In Fig. we show the largest eigenvalue of the 
static spin (or charge) susceptibility x'CQjO) m the Bril- 
louin zone, for the isotropic case t'/t = 1.0. For the 
original Kagome lattice, the susceptibility should have 
the exact six- fold rotational symmetry in the BZ. How- 
ever, since in our treatment the original Kagome lattice is 
transformed into the topologically-equivalent decorated 
square lattice (isotropic Kagome lattice), the susceptibil- 
ity has a distorted symmetry. Then the largest values 
of the susceptibility are located at q « ±(0.47r, — OAtt), 
±(0.47r, 0.87r), and ±(0.87r, OAtt) and this is due to nest- 
ing behavior. Although there is not the strong nesting 
because of the rather isotropic Fermi surface, when the 
Fermi surface touches itself in a BZ at several points for 
a given nesting vector, this enhances the susceptibility. 
Such conditions are illustrated in Fig. 0] (a) , where the 
arrows denote the corresponding nesting vectors. It is 
seen that these vectors indeed give rise to large values of 
the susceptibility. The results for the anisotropic cases 
are also shown in the left panels of Fig. and they can 
also be explained by the Ferm surface nesting. 

In the following, we will study the effects of geometrical 
frustration of the isotropic/anisotropic Kagome lattice 
with including the Coulomb interaction U. Of course, 
the strength of electron correlation is characterized by 
the ratio of the Coulomb interaction to the kinetic energy, 
but the latter is insensitive to the hopping anisotropy t'/t 
and therefore the strength is essentially given by the ratio 
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FIG. 4: (a) Nesting of Fermi surface. Each shaded ellipse 
depicts a Fermi surface in the BZ shown by a square box. 
Arrows show the nesting vectors, (b) Contour plot of the 
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FIG. 5: The t dependence of the total bandwidth at non- 
interacting system. The inset shows the magnification. 

U /t. This is because the change in the total bandwidth 
and the DOS at the Fermi energy is small. We show the 
change in the bandwidth as a function of t' in Fig. 
Although the bandwidth decreases slightly with decreas- 
ing hopping t' , its change is very small. Even at t'/t = 0, 
the total bandwidth is 4\/2t as 5.66i, which is about 94% 
of that of the isotropic case (t'/t — 1.0). This enables us 
to systematically discuss electron correlations due to U 
by using the anisotropic model. 



B. Magnetic properties 

We now turn to study electron correlations in the 
isotropic/anisotropic Kagome lattice, and let us start 
with magnetic properties. Shown in Fig. is the 
Coulomb interaction dependence of the maximum eigen- 
values of the static spin susceptibility X^(q, 0) for three 
typical values of t' . We first discuss the isotropic case 
(t'/t = 1.0). For U/t = 0, the largest values of the 
susceptibility are located at the 6 points in the BZ, 
q « ±(0.4vr,-0.47r), ±(0.4tt, 0.8tt), and ±(0.8tt, 0.4tt). 
With increasing U, the susceptibility is enhanced in the 
whole BZ. However, its q-dependence is reduced in con- 
trast to the unfrustrated cases. In particular, it is found 
that the susceptibility for various q's besides the above- 
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FIG. 6: The maximum eigenvalues of the static spin susceptibility % (*J)0) f° r three typical U values. Upper, middle, and 
lower panels correspond to t'/t = 1.0, 0.8, and 0.6, respectively. The temperature T/t is 0.1. 



mentioned 6 points are also strongly enhanced with in- 
creasing U, so that the q-dependence of the spin suscep- 
tibility is weakened and spin fluctuations become more 
isotropic in the BZ. We also show the susceptibility for 
anisotropic cases t'/t = 0.8 and 0.6 in Fig. |SJ Al- 
though the (distorted) six-fold symmetry is lost due to 
hopping anisotropy, the spin susceptibility at U/t = 
still has large values at specific q-positions, which are 
also able to be expected by the nesting picture similar 
to t'/t = 1.0 case. When we switch on U, the peaks 
in the q-dependence at U/t = gradually grows and 
this enhancement is more remarkable for smaller t' . For 
example, for t'/t = 0.6, the peak positions in the inter- 
mediate U region are q ~ ±(0.57r, — 0.57r), which implies 
that the system prefers the spin configuration with the 
period of twice the original unit cell. 

We note here that just at t'/t = 0, a flat band is located 
at the Fermi level at half-filling. It is known that this 



gives rise to the so-called flat-band ferromagnetic ground 
states. The nature of the ferromagnetism in this case 
has been well studied, so that we do not give detailed 
discussions for the cases t'/t <C 1.0. 

In order to discuss the instability to spin ordered 
states, it is instructive to compare the FLEX results 
with those calculated by the random phase approxima- 
tion (RPA) . The inverse of the largest eigenvalue of the 
spin susceptibility \ W(q, 0) in the whole BZ is shown 
in Fig. \7\ (a) for these two approximations. When this 
value becomes zero, the corresponding spin order, i.e., 
spin density wave (SDW), appears. Both results of the 
FLEX approximation and the RPA show that the dom- 
inant spin susceptibility is suppressed with increasing t' ', 
indicating the suppression of spin order. Note that the 
RPA spin susceptibility diverges at a finite U/t for all t', 
and the divergence is even enhanced, as seen in the con- 
cave shape of the lines, due to the multi-band structure. 
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FIG. 7: (a) The inverse of the largest eigenvalue of the spin 
susceptibility as a function of the Coulomb interaction U for 
various t' . Solid (dashed) lines represent the results of the 
FLEX approximation (RPA). (b) Largest eigenvalues of the 
spin and charge susceptibilities calculated by the FLEX ap- 
proximation as a function of U /t. Solid (dashed) lines denote 
the charge (spin) sector, respectively. All of the data are nor- 
malized by the value at U/t = 0.0. 



On the other hand, the enhancement in the FLEX results 
is much weaker at large U ft and the spin susceptibility 
does not seem to diverge. The behavior of the FLEX 
results is more reliable in the intermediate U jt region, 
since mode-mode couplings are completely neglected in 
the RPA and the instability is overestimated. In par- 
ticular, the suppression is most remarkable around the 
isotropic point t' ~ t, and we may conclude that no spin 
order is realized for the isotropic Kagome lattice within 
the FLEX approximation. 

Summarizing the magnetic properties, the spin sus- 
ceptibility is enhanced in the presence of the Coulomb 
interaction in the whole BZ, but geometrical frustration 
of the Kagome lattice has a strong effect to reduce the 
q-dependence of the susceptibility and to suppress the 
instability to any type of the SDW order. Although the 
present calculations are based on a weak coupling ap- 
proach, the obtained results are quite consistent with 
properties expected in the strong coupling regime, the 
frustrated Heisenberg spin system, where the magnetic 
order is suppressed due to strong frustration. 

We have so far discussed the instability of spin-triplet 
polarization, but the instability in the spin-singlet chan- 
nel is even weaker. The largest eigenvalues of the charge 
susceptibility x^(q, 0) m the BZ are always smaller than 



that for %( s ) (q, 0) within the FLEX approximation and 
there is no indication of charge order (charge density 
wave, CDW) or current-carrying state (shown in Fig. 
(b)). Thus we do not give detailed discussions on the 
charge susceptibility here. 

C. One-Particle Spectral Function 

We now discuss the one-particle spectral function and 
the total spectral density, which are defined by 

A(k,u) = --!- Y]lmG mm (k,w + i£) (13) 

m 

= ^E^ k ' w )' ( 14 ) 

k 

where 6 is a small adiabatic constant. 

In order to obtain dynamical quantities in the FLEX 
approximation, the analytic continuation from the imagi- 
nary Matsubara frequency to the real frequency is usually 
performed by using the Pade approximation or the max- 
imum entropy method. However, it is known that this 
procedure sometimes encounters numerical difficulties. 
In the present study, we directly calculate the spectral 
function in a real frequency formulation without resort- 
ing to analytic continuation. (Our formulation is similar 
to those used in the previous worksS^&Si, but simpler 
from the viewpoint of numerical techniques.) We need a 
small but finite adiabatic constant S for the Green func- 
tion, but its effects can be corrected systematically and 
we have checked that the Pade approximation^ repro- 
duces similar results. We show here the results calculated 
by our real frequency technique, since they are free from 
numerical errors due to analytic continuation. 

Figure [3] shows the one-particle spectral function for 
several k's near the Fermi surface. All the spectra shown 
consist of three peaks reflecting the three-band structure; 
namely the cv-th peak is mainly due to the contribution 
of the band with E^ a . As U increases, two peaks at 
large \u\ is gradually smeared due to life-time effects. 
On the other hand, the peak near u> = keeps a sharp 
structure, implying a well-defined Fermi liquid behavior. 
This behavior is particularly prominent in the isotropic 
case t'/t = 1.0 irrespective of k's, (a)-(c). However in 
the anisotropic cases of t'/t = 0.8 and 0.6, the central 
peak changes its shape from (a) k rss (0.757T, 0.757r) to (c) 
k w (0.42-7T, — 0.427r) where the amplitude of the spectra 
at to ~ is reduced with the decrease of t'/t. This behav- 
ior suggests that quasi-particles are stabilized in the re- 
gion of isotropic hopping t' ~ t. As discussed in the previ- 
ous subsection, increasing t' represents the enhancement 
of effects of geometrical frustration, and this leads to the 
trend that the wavevector dependences are strongly re- 
duced in the BZ. Since the spin fluctuations are local- 
ized in the real space and these amplitude are strongly 
reduced, various orders, such as the SDW, CDW, and 
currant-carrying states, are suppressed and the metallic 
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FIG. 8: One-particle spectra for various k. (a) k 
» (0.757r,0.757r), (b) k » (0.75tt,0), and (c) k « 
(0.427T, — 0.427r), respectively, (d) Fermi surface for several 
values of t' . Small circles represent the location of k cho- 
sen for (a)-(c). The energy is measured from the Fermi level 
{u> = 0). The Coulomb repulsion U/t — 2.5, the temperature 
T/t = 0.1 and adiabatic constant 8 = 0.05. The insets show 
the magnification of the low-frequency regime. 




FIG. 9: Total spectral density for various U. (a) t'/t = 1.0, 
(b) t'/t = 0.8 and (c) t'/t = 0.6. 



state is stabilized. Therefore the quasi-particle picture 
becomes fine description. This result indicates that ge- 
ometrical frustration on the Kagome lattice may be im- 
portant for the formation of quasi-particles. 

For reference, we also present in Fig. [U] the results 
of the total spectral density defined by Eq. JT3J. As 
U increases, the total spectral density is smeared in the 
whole energy region due to the life-time effect. However, 
in the weak-coupling regime treated here, the total spec- 
tral density near the Fermi level is rather insensitive to 
the change in anisotropic hopping t' in comparison with 
other physical quantities, such as the one-particle spec- 
tra, etc. Although the one-particle spectral weight with 
k for the k y = —k x direction close to the Fermi surface 
is strongly reduced at uj ~ with decreasing t'/t, the 
amplitude of total spectral density at uj ~ becomes 
largest at t'/t = 0.6 even when U is introduced. In the 
intermediate U region, the one-particle spectrum at k far 
from the Fermi surface mainly consists of the incoherent 
parts. Since the local spectral density is obtained by the 
k-summation in Eq. 1)14(1 ■ the contributions of these in- 
coherent parts far from Fermi surface become quite large 
even at p(0) with increasing U . Therefore, it is difficult to 
discuss about the coherence of quasi-particles only based 
on the total spectral density p(0), and we need more de- 
tailed information from A(k, to). 

Summarizing this subsection, as compared with the 
total spectral density, the one-particle spectrum gives 
us more detailed information on the nature of quasi- 
particles. The weight of the coherent part in the one- 
particle spectrum grows with increasing t', which means 



that the quasi-particle behavior is stabilized by geomet- 
rical frustration. 



D. Renormalization effects 

To investigate the properties of quasi-particles in more 
detail, we calculate the renormalization factor and the 
Fermi velocity of quasi-particles. According to the con- 
ventional Fermi liquid theory, the one-particle spectrum 
consists of the (5-function like coherent part and the broad 
incoherent part. As the Coulomb interaction increases, 
the quasi-particle mass is enhanced and the weight of the 
coherent part becomes small. For the single band case, 
the renormalization factor is given by the self energy as 



<9Re£(k,w) 



du 



j=0- 



(15) 



which corresponds to the weight of the coherent part. 
However, in the present case, the one-particle spectrum 
at a given k has contribution of three bands. When the 
coherent quasi-particle part is separated from the other 
parts, it is given by 



A coh (k,uj) 



l 



7T (W - £ k ) 2 + 7 2 



+ b, 



(16) 



where £k is the shift from the Fermi level, b is the correc- 
tion due to the other peaks which are located far from 
the Fermi level, and 7 represents the inverse of life-time 
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of the quasi-particle. These parameters can be deter- 
mined by numerical fitting of the data of the spectral 
function. However instead of performing this procedure, 
we here use an alternative approach to simply determine 
the renormalization factor. Since the middle band near 
the Fermi level is isolated sufficiently far from the other 
bands, we can introduce the following effective self-energy 
in the low-energy region. 



5(k,w)- 1 -G(k,w)" 



where 



m 

G(k,w) = ~^G mm (k,u>). 



(17) 

(18) 
(19) 



The approximate formula for the renormalization factor 
then reads, 



1 



<9Re£ ctt (k,w) 



dw 



j=0- 



(20) 



Although Eq. i|20[l is not exact, this formula works well 
as far as the low energy regime is concerned. We have 
checked for several k that the results obtained by using 
Eq. Ij20|l are in good agreement with the renormalization 
factor determined by the fitting Eq. H16fl . 

Figure ITUl shows the k-dependence of the renormaliza- 
tion factor as a function of U for various t'. We show the 
results only in the k region, < \k y \ < k x <ir, because 
the values in other regions are readily obtained accord- 
ing to the symmetry property. In all the cases of t'/t, 
the renormalization factor decreases with the increase 
of U, meaning that the effective mass is enhanced. An 
important point is that becomes nearly independent 
of k in the isotropic case (t'/t = 1.0), being consistent 
with the results of one-particle spectral functions. On 
the other hand, with decreasing t' , the difference in Z^ 
for different k's becomes more remarkable. In particular, 




FIG. 10: Renormalization factor as a function of U at tem- 
perature T/t = 0.1; (a) t'/t = 1.0, (b) t'/t = 0.8, and (c) 
t'/t = 0.6, respectively. The insets show the Fermi surface 
(solid line) in the half BZ (0 < k x < ir) where symbols (open 
circles etc.) denote various k for which is calculated. 



the Zk for k's along the direction of k y — —k x is strongly 
renormalized. It should be noted that the largest Z-^ in 
the anisotropic case (open squares in Fig. EII(k) and (c)) 
has almost the same amplitude of that for t'/t = 1.0. In 
other words, with decreasing t' , the quasi-particle weight 
is reduced in comparison with that of t'/t — 1.0. There- 
fore this is another evidence that geometrical frustration 
stabilizes the metallic state up to larger U regime when 
t'/t = 1.0 (the isotropic Kagome lattice). 

Shown in Fig. ^2 is the absolute value of the Fermi 
velocity, which is given by 



dk x 



dk v 



(21) 



which is obtained from the renormalized quasi-particle 
energy of the middle band, 



Eu = Z\, 



Re£ ctt (k,0) 



(22) 



In the case of t'/t ~ 1.0, the Fermi velocity has the dis- 
torted six-fold rotational symmetry in the BZ reflecting 
the lattice structure for (Fig. El(d)) and the ratio of the 
largest and smallest velocity is about 1.17 for U/t = 0. 
With increasing U, the ratio is gradually reduced to the 
smaller value of about 1.10 for U/t = 2.0, indicating 
again the suppression of the k-dependence. On the other 
hand, for the anisotropic case t'/t < 1.0, the ratio is 
not reduced but enhanced. For instance, for t'/t = 0.6, 
the ratio changes from about 1.79 for U/t = 0.0 to 2.45 
for U/t — 2.0. These results are in accordance with 
those of the renormalization factor in the sense that the 
wavevector dependence grows with increasing hopping 
anisotropy. 



(d) U/t =0.0, f/f =1.0 



> i 



o 



o 

-m 







U/t =0.0 


(a) t'/t =1.0 


1.0- 




2.0 










(b)f'A=0.8 

I.I.I 


1 1 1 1 1 1 1 

. (c)f'/f=0.6 /S'~ 











1.95 



1.80 



1.65 




(.T..T) 



K (7l,-7t) 

(e) 



/ o 









FIG. 11: Angular dependence of the Fermi velocity for various 
U; (a) t'/t = 1.0, (b) t'/t = 0.8, and (c) t'/t = 0.6 (Eq. lETlh 

(d) Fermi velocity at U/t = and t'/t = 1.0 in the k space. 

(e) angle cj> of the k point on the Fermi surface. 
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IV. SUMMARY AND DISCUSSIONS 

In this paper, we have studied how geometrical frus- 
tration of the isotropic/anisotropic Kagome lattice af- 
fects physical properties in a metallic phase. By using 
the FLEX approximation for the Hubbard model at half- 
filling, we have calculated the spin and charge suscepti- 
bility, the one-particle spectrum, the total spectral den- 
sity, the quasi-particle renormalization factor, and the 
Fermi velocity. In the isotropic Kagome lattice, we have 
shown that the spin susceptibility is enhanced by the 
Coulomb interaction as far as the amplitude of suscep- 
tibility is concerned. However, in contrast to ordinary 
cases of magnetic instability, its wavevector dependence 
is strongly suppressed at the same time, and there ap- 
pear many peaks in the BZ accumulating to form a line. 
This reduction of the wavevector dependence around the 
Fermi surface is also observed for other quantities, such 
as the quasi-particle weight of the one-particle spectrum. 
We have demonstrated that this behavior is indeed due to 
the frustrated lattice geometry by controlling frustration 
by anisotropic hopping; the k-dependence is considerably 
recovered with the decrease of t' i.e. when frustration is 
reduced. By comparing the spin susceptibility calculated 
by the FLEX approximation with that of RPA, we have 
shown that the instability to magnetically ordered phases 
is dramatically suppressed by mode-mode couplings of 
fluctuations, which are neglected in the RPA level. This 
suppression is most prominent in the isotropic case, i.e. 
most frustrated case. Furthermore, considering the one- 
particle spectrum, we have shown that the coherent part 
has a large weight in the isotropic Kagome lattice up to 
larger U regime, which is also consistent with the be- 
haviors of the renormalization factor and Fermi velocity. 
These results indicate that geometrical frustration stabi- 
lizes the formation of quasi-particles. 

The main conclusion of the present study is that the 
strong geometrical frustration reduces the wavevector de- 
pendence of spin fluctuations and suppresses the insta- 
bilities to various orders, such as the SDW, CDW, and 
current-carrying states, so that the metallic state is sta- 



bilized and the quasi-particles become well-defined in the 
isotropic Kagome lattice. This conclusion derived from 
the FLEX approximation, which is a essentially weak- 
coupling approach, but it is consistent with the results 
known for the strong coupling model, i.e., the absence 
of any magnetic order due to strong frustration in the 
antiferromagnetic spin system on the Kagome lattice. 
Therefore, it is reasonable that characteristic properties 
discussed in the present work for the weak-coupling Hub- 
bard model on the Kagome lattice are naturally con- 
nected to those inherent in the frustrated Heisenberg spin 
system. 

Within the FLEX approximation in this paper, we 
have not observed other ordered states, either such as 
the staggered flux (current-carrying) state 4 * 4 *^*^ and the 
charge ordered stated etc., which are other possible or- 
ders in the strong coupling regime. In the present case of 
half-filling, although charge ordered states are unlikely, a 
staggered flux state may appear due to strong geometri- 
cal frustration. Furthermore, upon finite hole doping, ge- 
ometrical frustration may realize a charge ordered state. 
It is an interesting open problem to be addressed in the 
future. 

Experimentally, most compounds studied so far have 
been restricted to spin systems on the Kagome lattice, 
such as SrCrsGa^ m 48 ' 49 . It is thus interesting to study 
electron correlations in metallic Kagome lattice materi- 
als. 
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